Effect of structural defects on anomalous ultrasound propagation in solids during 

second-order phase transitions 
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The effect of structural defects on the critical ultrasound attenuation and ultrasound velocity 
dispersion in Ising-like three-dimensional systems is studied. A field-theoretical description of the 
dynamic effects of acoustic-wave propagation in solids during phase transitions is performed with 
allowance for both fluctuation and relaxation attenuation mechanisms. The temperature and fre- 
quency dependences of the scaling functions of the attenuation coefficient and the ultrasound velocity 
dispersion are calculated in a two-loop approximation for pure and structurally disordered systems, 
and their asymptotic behavior in hydrodynamic and critical regions is separated. As compared to 
a pure system, the presence of structural defects in it is shown to cause a stronger increase in the 
sound attenuation coefficient and the sound velocity dispersion even in the hydrodynamic region as 
the critical temperature is reached. As compared to pure analogs, structurally disordered systems 
should exhibit stronger temperature and frequency dependences of the acoustic characteristics in 
the critical region. 

PACS numbers: 64.60.Ak, 64.60.Fr, 64.60.Cn, 71.23.-k, 43.35.+d 



I. INTRODUCTION 

The progress in understanding the nature of critical 
phenomena is mainly related to the theoretical and ex- 
perimental studies of critical dynamics in condensed mat- 
ter. However, the descriptions of the nonequilibrium be- 
havior of systems during phase transitions still contain 
a number of unsolved problems. This is due to the fact 
that studying the dynamic properties of critical fluctua- 
tions, which have anomalously high amplitudes and slow 
damping, encounters problems that are more complex 
than the problems that arise when equilibrium proper- 
ties are described. Qualitatively, this is caused by the 
necessity of taking into account the interactions between 
order-parameter fluctuations and other longlived excita- 
tions. 

The dynamics of phase transitions contains a number 
of physically important processes that are determined 
by the behavior of a multispin correlation function and, 
thus, are particularly complex for a theoretical descrip- 
tion. These are, for example, thermal processes near the 
critical point in a liquidgas system, the attenuation of 
electromagnetic-field energy that accompanies magnetic 
resonance phenomena, and the anomalous attenuation 
and scattering of acoustic waves in various media during 
phase transitions. The latter processes are important, 
since they underlie resonance and ultrasonic methods of 
studying critical dynamics. 

The unique feature of ultrasonic methods is the fact 
that, at temperatures that are close to a second-order 
phase transition temperature in magnetic systems and 
systems with structural phase transitions, researchers 
detect both anomalously strong ultrasound attenuation 
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and an anomalous change in the ultrasound velocity, 
which can easily be observed in experiment (Figs.[]Ji,[T]D) 
0, Q ■ These phenomena are caused by the interaction of 
low-frequency acoustic oscillations with longlived order- 
parameter fluctuations, which produce a random force 
that disturbs normal acoustic regimes by means of mag- 
netostrictive spinphonon interaction. In this process, re- 
laxation and fluctuation attenuation mechanisms can be 
distinguished. The relaxation mechanism, which is due 
to a linear dynamic relationship between sound waves 
and an order parameter Q, manifests itself only in an 
ordered phase, where the statistical average of the order 
parameter is nonzero. Since the relaxation of the order 
parameter near a phase-transition point proceeds slowly, 
this mechanism plays an important role in the dissipation 
of low- frequency acoustic oscillations. The fluctuation at- 
tenuation mechanism, which is determined by a quadratic 
relation between the deformation variables in the Hamil- 
tonian of a system with order-parameter fluctuations, 
manifests itself over the entire critical-temperature range. 
To date, there exist a considerable number of works 
that deal with a theoretical description of the ultrasonic 
anomalies that appear in condensed matter during phase 
transitions 0, [H, [6| and give an adequate explanation of 
experimental results 0, S, Q • 

One of the most interesting and important problems 
from both experimental and theoretical viewpoints is the 
study of the influence of structural defects on the ultra- 
sound propagation characteristics in materials undergo- 
ing phase transformations. The structural disorder in- 
duced by impurities or other structural defects plays a 
key role in the behavior of real materials and physical 
systems. The matter of particular interest is the effect 
of frozen structural defects, whose presence can manifest 
itself in a random perturbation of a local transition tem- 
perature, as it occurs, for instance, in ferro- and antifer- 
romagnetic systems in the absence of an external mag- 
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netic field. The statistical features inherent in systems 
with a frozen disorder create considerable difficulties for 
both the analytical description of the behavior of such 
systems and the experimental methods of their investi- 
gation. According to the heuristic Harris criterion [lol ]. 
the effect of frozen point defects becomes noticeable and 
induces a new type of critical behavior if the critical ex- 
ponent for the heat capacity of a pure system is positive. 
As was demonstrated in previous studies, this criterion 
is only met for Ising-like systems. Therefore, the effect 
of point structural defects on the critical behavior is neg- 
ligible in systems with a multicomponent order param- 
eter, such as the XY model and the Heisenberg model. 
Therefore, one of the most challenging problems from a 
physical point of view is the study of the effect of struc- 
tural defects on the critical behavior of systems with a 
single-component order parameter, in which the presence 
of a structural disorder leads to a substantial change in 
the critical-behavior characteristics. 

However, the problem of the effect of structural de- 
fects on the characteristics of ultrasound propagation in 
materials that undergo phase transformations still re- 
mains unsolved because of the complexity of the theo- 
retical description of the four-spin correlations of order- 
parameter fluctuations, which determine acoustic charac- 
teristics. Pawlak and Fechner attempted to describe 
the effect of point defects on the ultrasound propagation 
parameters near a critical temperature using the first- 
order e expansion. However, as was shown in our work 
some mistakes crept in the description of this phe- 
nomenon in [Tll | ; in particular, they used wrong diagrams 
for taking into account the dynamic effects of the interac- 
tion of order-parameter fluctuations via a defect-induced 



field and they did not use diagrams that give a notice- 
able contribution to the attenuation coefficient. More- 
over, earlier investigations [ID, [3] based on the field- 
theoretical description of pure and disordered systems in 
two-loop and higher approximations with the use of the 
method of summation of asymptotic series demonstrated 
that the results obtained in the lowest-order e expansion 
can only be considered as a crude estimate, especially 
for disordered systems. Thus, the results obtained in 
require reevaluation using a more precise approach. 
For this purpose, in this work we performed a correct 
field-theoretical description of the effect of structural de- 
fects on the anomalous critical ultrasound attenuation 
and the anomalous change in the ultrasound velocity in 
three-dimensional Ising-like compressible systems with 
allowance for both the fluctuation [12] and relaxation 
attenuation mechanisms without using the e-expansion 
method. 



II. THE MODEL 

For phase transitions in compressible systems, the rela- 
tion between an order parameter and elastic deformations 
is an important factor. As was first shown in [15J, the 
critical behavior of compressible systems with a quadratic 
striction is unstable in regard to the relation between an 
order parameter and acoustic modes, and a first-order 
phase transition that is close to a second-order phase 
transition is realized. However, as was clarified in 
the conclusions made in [l5| are valid only at low pres- 
sures and, beginning with a certain threshold pressure, 
the deformation effects induced by an external pressure 
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change the order of the phase transition. 

The Hamiltonian of a disordered compressible Ising 
model can be written as 



H = H K 



H n 



Hi,, 



Hi, 



(1) 



The contribution of the deformation degrees of freedom 
is determined as 



H„, 
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where u a p(x) are the components of the deformation 
tensor and C£ arc the elastic constants. The use of 
an isotropy approximation for H e i is caused by the 
fact that, in the critical region, the system behavior 
parameters are determined by an isotropic fixed point 
of renormalization- group transformations, while the 
anisotropy effects are negligible [17| . The magnetic com- 
ponent H op is represented in the form of the Ginzburg- 
Landau Wilson Hamiltonian 



H op = I d d x 



ir 5 2 + i(V5) 2 + ^ 4 



(3) 



where S(x) is the spin order parameter, u is a posi- 
tive interaction constant, and To = (T — Tq c )/Tq c is the 
reduced phase-transition temperature. The Hi n t compo- 
nent determines spin-phonon interaction, 



Hint = I d 



fjo 



(4) 



where go is the quadratic-striction parameter. The effect 
of defects is taken into account by the term 



H lmp = J d d x [At(x)S 2 ]+ J d d 
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. (5) 



where random and Gaussian-distributed variables At(x) 
and h(x) determine local phase-transition temperature 
fluctuations and random stress fields, respectively. 

To perform calculations, it is convenient to use the 
Fourier components of the deformation variables in the 
form 



where e\(q) is the polarization vector. 

We then perform integration with respect to the off- 
diagonal components of the uniform part of the defor- 
mation tensor u£l, in the statistical sum (they are not 
essential for the critical behavior of the system in an elas- 
tically isotropic medium) and obtain a Hamiltonian for 
the system in the form of a functional for the spin or- 
der parameter S (q) and the normal coordinates of the 
deformation variables Q\ (q) 

H = X - [ d d q (t q + q 2 ) S q S- q + [ d d q q h q Q q , x + 
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The relaxation critical dynamics of compressible sys- 
tems is described by dynamic equations of the type of 
generalized Langevin equations, 



dH 
dS~ 



•£ q + r hs, 



Q q ,\ 



dH 
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where To and Do are bare kinetic coefficients; £ g (x, t) and 
r] q (x, t) are Gaussian-distributed quantities that have the 
character of a random force; and hs and hQ are the fields 
thermodynamically conjugated to the spin and deforma- 
tion variables, respectively. 

When solving the set of nonlinear equations ([8|) with 
Hamiltonian H (S, Q) iteratively, we can single out the 
elastic- variable response function D(q,tu), which is de- 
termined as 

D ( Q) w ) = ffip^l = [(Q 9>u>x Q- q ,- u , x )] , (9) 
ohn 



U a /3 
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+ V 1/2 £ u a p(q) exp (iqx) , (6) and the spin-variable response function G(q, uj) 



where q is the wavevector, V is the volume, li^l 
is the uniform deformation tensor, and u a p(q) = 
i/2 [q a up + qpu a ]. We introduce an expansion in terms 
of the normal coordinates, 



G(q,u) 



6[(S q ,*)] 
Sh s 



[(S qtU S- qt - u )] , (10) 



where (...) stands for statistic averaging over random 
Langevin forces, [...] stands for the averaging over the 
fluctuations of random fields Ar_ g and h q that are spec- 
ified by structural defects, and to is the characteristic 
ultrasonic vibration frequency. 



Using the Dyson representation, we present the G(q, lo) 
and D(q,u>) response functions in the form 

G~ x {q,u>) = Gq 1 (q,uj) + IL(q,(j), (11) 
D- 1 (q,cj) = D 1 (q,u;) + j:(q,cj). (12) 

The bare Go(q,u>) and Do(q,u>) response functions are 
determined as: 

Do(q,u) = (to 2 - a q 2 - iuiD a q 2 ) 1 , 
G a (q,uj) = (ico/T +(T + q 2 )y 1 . 

In the low-temperature phase, the response function 
contains an additional relaxation contribution 



Sq = MSgfl 



with the magnetization 
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T <T. 



(13) 



(14) 



where B is the phenomenological relaxation parameter 
and tfq is the fluctuation part of the order parameter. 

The self-energy part E(<?, lo) of the D(q,Lo) response 
function is directly related to the dynamic characteristics 
of ultrasound propagation [l8| . 

As a result, the ultrasonic attenuation coefficient can 
be expressed through the imaginary part of £(g,w) 



a(u>, t) ~ lo ImS(0, lo), 



(15) 



and the sound velocity dispersion is expressed through 
its real part, 



2 (lo, t) - c 2 (0, t) ~ Re (S(0, lo) - E(0, 0)) . 



(16) 



We calculated T,(q,Lo) in a two- loop approximation. The 
diagrammatic representation of £(?, lo) is shown in Fig. [5] 
These Feynman diagrams contain (i-dimcnsional integra- 
tion (in our case, d = 3). 

When approaching the critical point, the correlation 
length £ tends to infinity, and, when £ _1 <C A (where A is 
the cutoff parameter of the integration over wavevectors) , 
the system characteristics demonstrate their asymptotic 
scaling behavior for wavevectors g« A. Thus, the calcu- 
lation of these quantities can be carried out in the limit 
A — * oo. The application of a renormalization-group pro- 
cedure eliminates the divergences that appear in the ther- 
modynamic variables and kinetic coefficients at A — > oo. 

To calculate attenuation coefficient (fT5")) and ultra- 
sound velocity dispersion (|16p and to eliminate the diver- 
gences in H(q, lo) at q — > we used the matching method 
in 19], which was then generalized for the description of 
the dynamic behavior of a system in [20] . 

Thus, using a scaling relationship for the dynamic re- 
sponse function 

D (q, lo, r) = e^ l D (qe l , (lv/T ) e* 1 , re l ' v ) , (17) 
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FIG. 2: Diagrammatic representation of S(g, lo) in 
a two-loop approximation. Solid line (a) corresponds 
to Go(q,Lo); lines (b) with a cross, to Ca(q,io) = 

2D7 1 ^/r ) 2 + (q 2 + r ) 2 j ; vertex with dashed line (c), 

to v = [(At) 2 ] ; and line (d), to relaxation insertion M 2 8 q fi. 



we can calculate the right-hand side of the equation for 
some constant value I = /*, at which not all of the argu- 
ments in the response function disappear simultaneously. 
The choice of I* is determined by the condition 



(a /T ) e 



zl* 



i/z 



+ q 2 e 21 ' 



1. 



(18) 



malization- group transformation, namely, to find a rela- 
tion between the behavior of the system in the precritical 
regime at a low value of reduced temperature r and the 
behavior of the system in a regime far from the criti- 
cal mode, i.e., without divergences in £(<;, lo). As was 
demonstrated in [201 ] . matching condition (|18p provides 
an infrared cutoff for all diverging quantities. Based on 
Eq. (fTS]) . we find the solution for I* in the form of a func- 
tional dependence on lo and r, which is specified by the 
static critical exponent v of the correlation length and by 
the dynamic critical exponent z, 



1 + (y/2) 



4/z 



-1/4 



(19) 



where y = lot~ zv /To is the argument of the F(y) func- 
tion. 

As is known from the theory of ultrasound scatter- 
ing in solids near a phase-transition temperature [H, [|| , 
the expression for the imaginary part of r) in the 
asymptotic limit (r — > 0, lo — > 0) can be defined by a 
scaling function <f>{y) 



Im£(u;, r) /lo 



(20) 



which depends on the single generalized variable y. At 
the same time, for the imaginary component of the self- 
energy part, the following scaling relationship is valid (l8j 



ImE(u) 



J(a+zv)/v 



lmS>(uje zl ) 
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The substitution of e l from Eq. (fT9)) into the right-hand 
side of this expression allows calculating the <f>(y) scaling 
function. 

In the asymptotic limit r — > 0, u — > 0), the expres- 
sion for the real part of E(o>, r) can be determined using 
another scaling function f(y), 

Re (E(0, w) - E(0, 0)) = r"" (/(t/) - /(0)) . (21) 

The real component of the self-energy part satisfies the 
scaling relation 

Re(E(0,£j) - E(0,0)) = 

= e' Q /"Re(E(0,u;e 2 ')-E(0,0)). (22) 

The dynamic scaling functions calculated in the two- 
loop approximation has the form 



2 r F a / u+l / 2v ~ 

7T y 2 



1 - 
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where g* , u* and t>* are the magnitudes of the interaction 
vertices at the fixed point of the renormalization-group 
transformations that corresponds to the critical behav- 
ior of the disordered compressible Ising model [2lj |. The 
terms in Eqs. (f2"3"]) and (pM]) that are proportional to M 2 
describe the relaxation contribution for the scaling func- 
tions of the attenuation coefficient and the sound velocity 
dispersion. In our subsequent numerical calculations of 
the scaling functions, we used the value v — 0.70 from [2lJ 
for the corresponding fixed point. The value of the dy- 
namic exponent [z = 2.1653) was taken from [22j], where 
the critical dynamics of a disordered Ising model was an- 
alyzed within the framework of a relaxation model. The 



use of this value of exponent z is valid in the case of dis- 
ordered Ising-like systems with a negative heat-capacity 
exponent, since the relation between an order parame- 
ter and elastic deformations in the critical dynamics of 
a compressible system exerts no substantial influence on 
the relaxation properties of the order parameter. 



III. ANALYSIS OF RESULTS AND 
CONCLUSIONS 

Perturbation-theory series are known to be asymp- 
totic, and the vertices of the interaction of order- 
parameter fluctuations in the fluctuation range r — > 
are too high to provide the direct application of Eqs. (|2"3"|) 
and ([2"4"|) . Therefore, to extract the necessary physi- 
cal information from the derived expressions, we apply 
the PadeBorel method, which is used to sum up asymp- 
totic series, that was generalized to a three-parameter 
case. Then, the forward and inverse Borel transforma- 
tions have the form 

oo 

<fi(w,u,v) = CijkW 1 u? v k = J e^ t B(wt,ut,vt)dt, 

i,j,k 

B(w,u,v) = J2 T. , lJk L M , w l u j v k , 
i,j,k {i+J+k)\ 

(25) 

where w = g 2 . 

To analytically continue the Borel transform of the 
function, we introduce a series in an auxiliary variable 
A 



B(w, u, v, A) 



E 



k k—i 
i=0 j=0 



, , — -w u 3 v \ 
k\ 



(26) 

and substitute it to the Pade [L/M] approximation at the 
point A = 1. This procedure was proposed and approved 
in [23| to describe the critical behavior of a number of 
systems containing several vertices of the interaction of 
order-parameter fluctuations. The fact [23| that the sys- 
tem retains its symmetry during the application of the 
Pade approximants in variable A becomes substantial for 
the description of multivertex models. In this work, we 
calculated the scaling functions in the two-loop approxi- 
mation using approximant [1/1]. 

The behavior of the dynamic 4>{y) and f(y) scaling 
functions calculated with summation methods for pure 
and disordered systems is shown in Figs. [3Ja) and [3£b) 
on a loglog scale. Depending on the interval of chang- 
ing variable y , the following asymptotic regions can be 
distinguished in the behavior of (f>(y) and f(y) : a hy- 
drodynamic region, where y ~ uj^ z ~ (q(,) z <C 1, and 
a critical region y ~ lu^ z 3> 1, which determines the 
behavior of the system near the phasetransition temper- 
ature (t = (T — T c )/T c -C 1). As is seen from these 
curves at y <C 1, the presence of a structural disorder 
does not affect the behavior of the <p{y) and f(y), scaling 
functions and, consequently, the behavior of this system; 
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however, it begins to manifest itself in the crossover re- 
gion 10 _1 < y < 10 and exerts an essential effect in the 
critical region y > 10. 

As follows from Eqs. (fl~5|) and ([20]) , the attenuation 
coefficient can be expressed as 

a(w,T)~J t T- a -'"il>(v), (27) 

and, using Eqs. (fTo| and (f2"Tj) . we can write the relation 
for the sound velocity dispersion in the form 

cV,T)-c 2 (0,T)~T-«(/(y)-/(0)). (28) 

The results of the calculations of the asymptotic de- 
pendences of the attenuation coefficient and the sound 
velocity dispersion for the critical and hydrodynamic re- 
gions are given in the table. The characteristics of their 
frequency and temperature dependences were determined 
in the range 10 -3 < y < 10 -1 , for the hydrodynamic 
regime and in the range 10 < y < 10 3 . for the critical 
regime. Note that, according to~[2p| . the real temper- 
ature range 10~ 3 < r < 10 _1 in ultrasonic studies of 
phase transitions corresponds to the range 1 < y < 10 2 , 
i.e., it covers the crossover region and the beginning of 
the critical region (precritical regime). 

It follows from the table that anomalously strong ul- 
trasound attenuation should be observed in both pure 
and disordered systems. For the disordered model, 
the increase in the attenuation coefficient as the phase- 
transition temperature is approached is expected to be 
stronger than that in the pure model even in the hy- 
drodynamic region, whereas in the critical region, the 
disordered system should exhibit stronger frequency and 
temperature dependences of the attenuation coefficient 
as compared to the pure system. 

These conclusions are supported by the model repre- 
sentation of the results of the numerical calculations of 



the critical temperature behavior of the attenuation co- 
efficient for both the pure and disordered systems per- 
formed at B = 0.3 and lo/Tq = 0.0015. (Fig. [J). These 
values were determined when we compared the calcu- 
lated temperature dependence of the attenuation coeffi- 
cient and the results of experimental studies of pure FeF2 
samples (Fig. |4l (dots) [ij, which demonstrate Ising-like 
behavior in the critical region. 

An analysis of the data related to the sound velocity 
dispersion (see table) demonstrates that, as compared to 
pure analogs, a structural disorder in Ising-like systems 
leads a stronger temperature dependence of the sound 
velocity dispersion in both the hydrodynamic and crit- 
ical regions and is characterized by an increase in the 
exponent of the temperature dependence (a decrease in 
the absolute value) when going from the hydrodynamic 




' 

Tt,G 72,0 72,4 72,8 73,2 72-.S 77,6 78.0 78,4 78,8 

T, K 

FIG. 4: Temperature dependence of the attenuation coeffi- 
cient calculated for (1) pure and (2) disordered systems at 
B — 0.3 and to /To = 0.0015. (dots) The results of experimen- 
tal studies of pure FeF2 samples [l[. 
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TABLE I: Asymptotic behavior of the sound attenuation coefficient and the sound velocity dispersion in the critical, precritical, 
and hydrodynamic regimes for pure and disordered systems 













Attenuation coefficient 










Regim 








Pure 




Disordered 












T <T C 


T>T C 


T<T C 


T>T C 


Critical 
Precritical 
Hydrodynamic 




y = 
y = 
y = 


10 1 - 
10 1 - 
10 -3 


10 3 
10 2 

-lo- 1 


w 0.98 T -0.08 

w i-08 T -o.ai 


W 1.05 T -0.17 
w 1.2Q T -0.37 


W 1.12 T -0.10 
w 1.22 r -0.25 

W 2 T" L44 


w 1.21 T -0.24 
w 1.37 r -0.48 


Sound velocity dispersion c 2 (lj,t) — c 2 (0, r) 




Regim 








Pure 




Disordered 












T < T c 


T>T C 


T < T c 


T>T C 


Critical 
Precritical 
Hydrodynamic 




y = 
y = 
y = 


10 1 - 
10 1 - 
10 -3 


10 3 
10 2 

-nr 1 


w o.u T -o.a6 

w 0.30 r -0.49 


w 0.34 r -0.54 
w 1.08 r -1.48 


W 0.26 T -0.31 
w 0.41 r -0.54 

C^T" 2 - 95 


w 0.49 r -0.66 
w 1.01 r -1.45 



to the critical region. However, the effect of a structural 
disorder has the converse character for the exponent of 
the frequency dependence of the sound velocity disper- 
sion: in the hydrodynamic region, the exponents of these 
two types of systems coincide, whereas, in the critical 
region, the sound velocity dispersion of a disordered sys- 
tem has a stronger frequency dependence compared to 
the pure system. The exponent decreases strongly when 
going from the hydrodynamic to the critical region. 

A particularly important result of our investigation 
consists in the predicted manifestation of the dynamic 
effects of structural defects on anomalous sound attenu- 



ation and sound velocity dispersion over a wider temper- 
ature range near the critical temperature (already in the 
hydrodynamic region) in comparison with other experi- 
mental methods [2J] , which require a narrow temperature 
range (of about r ~ 10~ 4 . to be studied for revealing 
these effects. Thus, the results obtained can serve as 
a reference for purposeful experimental investigations of 
the dynamic effects of structural defects on the critical 
behavior of solids using acoustic methods via the detec- 
tion of the influence of structural defects on the frequency 
and temperature dependences of the ultrasound attenu- 
ation coefficient and ultrasound velocity dispersion. 



[1] A. Ikushima, R. Feigelson, J. Phys. Chem. Solids. 32, 
417 (1971). 

[2] Kh. K. Aliev, I. Kh. Kamilov, and A. M. Omarov, Zh. 
Eksp. Teor. Fiz. 95, 1896 (1989) [Sov. Phys. JETP 68, 
1096 (1989)]. 

[3] L. D. Landau and I. M. Khalatnikov, Dokl. Akad. Nauk 

SSSR 96, 496 (1954). 
[4] A. Pawlak, Phys. Rev. B 44, 5296 (1991). 
[5] A. M. Schorgg and F. Schwabl, Phys. Rev. B 49, 11682 

(1993). 

[6] I. K. Kamilov and Kh. K. Aliev, Usp. Fiz. Nauk 168, 953 

(1998) [Phys. Usp. 41, 865 (1998)]. 
[7] R. A. Ferrel, B. Mirhashem, and B. Bhattacharjee, Phys. 

Rev. B 35, 4662 (1987). 
[8] T. J. Moran and B. Liithi, Phys. Rev. B 4, 122 (1971). 
[9] M. Suzuki and T. Komatsubara, J. Phys. C 15, 4559 

(1982). 

[10] A. B. Harris, J. Phys. C 7, 1671 (1974). 

[11] A. Pawlak and B. Fechner, Phys. Rev. B 40, 9324 (1989). 

[12] P. V. Prudnikov, V. V. Prudnikov, J. Phys.: Condens. 

Matter. 17, L485 (2005). 
[13] V. V. Prudnikov, P. V. Prudnikov, A. A. Fedorenko, 

Phys. Rev. B 62, 8777 (2000); 63, 184201 (2001). 
[14] R. Folk, Yu. Golovach, and T. Yavorskioe, Usp. Fiz. Nauk 



173, 175 (2003) [Phys. Usp. 46, 169 (2003)]. 

[15] A. I. Larkin and S. A. Pikin, Zh. Eksp. Teor. Fiz. 56, 
1664 (1969) [Sov. Phys. JETP 29, 891 (1969)]. 

[16] Y. Imry, Phys. Rev. Lett. 33, 1304 (1974). 

[17] Yu. A. Izyumov and V. N. Syromyatnikov, Phase Tran- 
sitions and Crystal Symmetry (Nauka, Moscow, 1984; 
Kluwer, Dordrecht, 1990). 

[18] H. Iro, F. Schwabl, Solid State Commun. 46, 205 (1983). 

[19] D. R. Nelson, Phys. Rev. B 14, 1123 (1976). 

[20] R. Folk, H. Iro, F. Schwabl, Z. Phys. B 27, 169 (1977). 

[21] V. V. Prudnikov and S. V. Belim, Fiz. Tverd. Tela (St. 
Petersburg) 43, 1299 (2001) [Phys. Solid State 43, 1353 
(2001)]. 

[22] V. V. Prudnikov, S. V. Belim, A. V. Ivanov, et ah, Zh. 
Eksp. Teor. Fiz. 114, 972 (1988) [Sov. Phys. JETP 87, 
527 (1988)]. 

[23] K. B. Varnashev and A. I. Sokolov, Fiz. Tverd. Tela (St. 
Petersburg) 38, 3665 (1996) [Phys. Solid State 38, 1996 
(1996)]; A. I. Sokolov, K. B. Varnashev, and A. I. Mu- 
drov, Int. J. Mod. Phys. B 12, 1365 (1998); A. I. Sokolov 
and K. B. Varnashev, Phys. Rev. B 59, 8363 (1999). 

[24] N. Rosov, C. Hohenemser, and M. Eibschutz, Phys. Rev. 
B 46, 3452 (1992). 



